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Abstract—Mean filters are among the most commonly used 
filters in the field of digital signal processing in general and image 
processing in particular. Therefore, the research and improve- 
ment of the filter technology will have a positive impact on both 
the hardware and the software. Our studies in this paper have 
shown a new implementation solution, which greatly improves 
computational complexity, thus reducing the time taken for 
implementation to be superior to that of the original technique. 
in experimental results in the form of software. In particular, 
these suggestions can naturally expand into the n-dimensional 
space. 


Index Terms—2-D mean filter, Noise filter, Blur filter, Digital 
image processing, Digital signal processing 


I. INTRODUCTION 


Mean filter in two-dimensional (2-D) space, also known 
as “2-D Mean Filter,’ is a filter commonly used in the field 
of digital image processing to perform tasks of smoothing, 
blurring, sharpening details, or removing noise [2], [3], [41], 
[5] and many other applications in the field of digital signal 
processing. Therefore, research is needed to improve imple- 
mentation techniques of the filter [6], [7], [8], as the results 
will contribute to simplify circuit architecture of the digital 
signal processing (DSP) for hardware solutions, or improve 
the filtering execution time for software solutions. This paper 
focuses on existing implementation techniques of the mean 
filter, and propose improvements to bring about a better 
possible outcome than current ones. 

The effect of 2-D mean filtering on mandil_color image 
(Figure |[I{a)) with a 5 x 5 2-D mean filter is shown in Figure 
[I{b) and that with 9 x 9 is shown in Figure [I{c). 





Figure 1: 2-D mean filter: (a) Original mandil_color image; 
(b) 2-D mean filtering with filter size = 5 x 5; (c) 2-D mean 
filtering with filter size = 9 x 9. 


II. PRELIMINARIES 


This section explains some of the important material related 
to 2-D mean filtering. 


A. 2-D Convolutions and Separable Filters 


In the field of digital signal processing, a digital image X 
of size MI x N, where M and N represent number of columns 
and rows respectively, is considered a 2-D signal x(n1, 72), 
where: 


( ) X(n1,N2) fO<n <M—-1,0<n<N-I1, 
Tiss) = , 
le 0 otherwise, 


and the 2-D filter H of the size m x n is determined by its 
impulse responses h(n 1, N2). The 2-D filter H is called a finite 
impulse response (FIR) filters if and only if: 


A(ny,n2)>0 if0<n<m0<ng<n 
h(n,,n2) =O otherwise 


where 0 < m< M—1and0<n< N—1. When m = 7, 
the filter is said to be square. 

Filtering of the input image X with the filter H is im- 
plemented using a 2-D convolution, which is given by the 
expression 


) * h(n, n2) 
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The output is a digital image Y represented by 2-D signal 
y(n 1,2) given by 


y(ny,noq)>0 if0<n <M+m—-10<n<N4+n-1, 
y(n1,n2) =0_ otherwise. 


From this it follows that the size of the output image will be 
(M+m-—1)x(N+n-1), and (N+ M-—-1)x(N4n- 
1) x (m x nm) arithmetic operations (one arithmetic operation 
= one multiplication and one addition) are required. Since 
the size of the original image is significantly larger than the 
size of the filter G.e., M,N >> m,n), then approximately 
(N x M) x (n x m) arithmetic operations are required. 

If h(n,, 2) is a separable sequence, it can be expressed as 


h(n1, 2) = hy(n1)ho (no) (2) 


where 


h1(n1) =0 outside O< 
he(n2) = 0 <n 


outside 0 
From (I) and (2), 


y(ni,n2) = ee m2) * ae ng) 
= $3 ae (ky, ka) hy ( (n1 — ky )ho(ne — kg) 
ky =—0o kg=—00o 
= Ss hy(n, — ky) S x(k, k2)he(n2 — ke). 
ky =—0o kg=—0o 
(3) 
For a fixed kj, 
Ss” x(k, kz )ho(ne = kz) 
kg =—0o 
in is a 1-D convolution of x(k1,n2) and h2(n2). For 
example, using the notation 
f(ki,n2)= > 2(ki,k2)ho(n2 — ke), (4) 
kg=—0o 


f (0,2) is the result of 1-D convolution x(0, n2) with hg(n2), 
as shown in Figure |2| Since there are M different values of ky 
for which x(k,,k2) is nonzero, computing f(k , 72) requires 
M 1-D convolutions and therefore requires approximately 
(N+n-—1)xnx M arithmetic operations. Once f(k1, 2) 
is computed, y(nm1, 2) can be computed from (3) and (4) by 


y(n1, 22) 


sy hy(n1 — ky) f (ki, n2). (5) 


ky=-—co 


From (5), for a fixed ng, y(n1,N2) is a 1-D convolution of 
hy(n,) and f(n1, 2). For example, y(n 1, 1) is the result of a 
1-D convolution of f(m 1,1) and hy(n 1), as shown in Figure 
where f(n1,7n2) is obtained from f(k1,n2) by a simple 
change of variable. Since there are (N+n—1) different values 
of ng, computing y(n1, 2) from f (ki, 2) requires (N +n— 
1) 1-D convolution thus approximately (IV + m-—1) x m x 
(N +n-—1) arithmetic operations. Computing y(n1, 2) from 
and (5), exploiting the separability of h(n1,n2), requires 
(N+n—-1)xnx M]4+|[(M+m-1)xmx(N+n-1)] 
arithmetic operations. Since the size of the original image is 
significantly larger than the size of the filter (.e., M,N > 
m,n), then approximately (N x M) x (n x m) arithmetic 
operations are required. 


B. 2-D Mean Filter 


A 2-D mean filter of size m x n is defined by the impulse 
response given by the following formula: 
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Figure 2: Example of convolving x(n1, 2) with a separable 
sequence h(n1, 2). 





: fO<ny<m-1,0<n<n-Il, 


h(n1, 2) _ - 


0 otherwise. 
(6) 


Figure Bla) shows the values of impulse response h(n1, 2) 
of size 3 x 3 expressed as a matrix (also known as a filter 
mask), and its equivalent form containing the coefficients of 
all 1’s with a normalizing constant of 1/9 (also known as the 
Gain of the filter). This type of integer conversion provides the 
ability to implement a 2-D mean filter on integer processing 
systems. Since all coefficients of the filter mask are converted 
to values 1, computing y(n1,n2) as in does not require 
the multiplication x(k,,k2)h(n, — ky,ng — ko). Therefore 
implementing a 2-D mean filter on an input image signal by 
requires approximately (N x M) x (n x m) additions and 
(N x M) divisions. 


III. SPATIAL SEPARABLE CONVOLUTION (SSC) 
TECHNIQUE 


The 2-D mean filter may be realized in several ways. 
The basic technique is a direct-convolution realization using 
the relation (Ip. When the sizes of the image and the filter 
increase, the numbers of additions and divisions required by 
the basic technique increase significantly. It is interesting to 








(a) 3 x 3 mean filter. (b) 5 x 5 mean filter. 
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Figure 4: A 2-D mean filter mask can be decompsed into a 
horizontal and vertical filter mask: (a) 2-D mean filter mask 
Aynxn of the size m x n; (b) horizontal 1-D mean filter mask 
Ayn «13 (c) vertical 1-D mean filter mask H),.,. 
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note that a 2-D mean filter can be seperated into two 1-D 
mean filters according to the definitions at and (6). We 
call this the spatial separable convolution (SSC) technique. 
This technique is so named because it deals primarily with 
the spatial dimensions of an image and filter which are the 
width and the height, and it utilizes the separability of a 2-D 
mean filter. Specifically: 














1 1 
h(ni,n2) = ih on hi(ni)ha(n2) (7) 
where 
1 if0<ni<m 
hia) =< aera 8 
(71) . otherwise, - 
and 
1 if0<ng<n 
hg (tis pS Kee = * 9 
2(n2) f otherwise. ) 


An example of a 3 x 3 mean filter is shown in Figure [3[a). 

Thus, a 2-D mean filter of any size can be quickly im- 
plemented via two 1-D filters according to the procedure 
described in Figure |2| In this way, we infer that the number 
of mathematical operations to perform will decrease from 
(N x M) x (n x m) additions and (NV x M) divisions in basic 
technique to (NV x M) x (n+ m) additions and (N x M) 
divisions. In other words, when m = n, the number of 
additions is being reduced is (N x M x m)/2. 

Figure 4 shows the filter mask of a 2-D mean filter of the 
size m xX nm can be decomposed into a sequence of two 1- 
D filters in different directions (usually horizontal, and then 
vertical). 


IV. PROPOSED VERY FAST MEAN FILTERING (VFMF) 
TECHNIQUE 


Our very fast mean filtering (VFMF) technique improves 
upon the SSC technique by utilizing the fact that major addi- 
tion computations are repeated as the filter moves from point 
to point in an image. To begin, let us focus on analyzing some 
important characteristics in the implementation of convolving 
an input signal 7(n1,2) with two component filters h(n1) 
and h2(n2) which are decomposed from the separable filter 
h(n,,n2), as presented in Section 


To begin, we perform a 1-D convolution first for each 
column of the 2-D image signal x(n1,n2) with the impulse 
response h2(nz) to obtain the intermediate signal f(k,,n2). 
Then (4) becomes 
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x(ky, ke). (10) 


For every value 7 of the variable m2, we have 
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De 


kg=i1-n+1 


f (ki, 7%) = x(k1, ke), 


and so f(k1,i +1) can be computed recursively using 


atl 
kg =1-n+2 


a 


L(ky, ke) 


= S- a(ki,ke)+0(k1,i+1) —2(ky,i-n+1) 


kg=i-n+1 
f(ki, 1) ae x(ky, 4 + 1) = x(k, (2 + 1) — n). 
(11) 


Thus, for each fixed value kj, the output of f at (k1,i+1) is 
calculated by adding the output of f at (k,,2) and the input of 
x at (ky1,i+1), and subtracting the value of x at (ky, (2+1)— 
n). So, obtaining each output f requires 2 additions, except for 
the first f of the column corresponding to nz = O, n additions 
are required. Figure |5| illustrates the recursive expressions of 
f’s. 

To continue, we then perform a 1-D convolution for each 
row of the intermediate signal f(k1,n2) with the impulse 
response hi(n1) to obtain the output y(n 1,2). Then 
becomes 
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Figure 5: Computation of f(k,,7+ 1) from f(k,,7). 


y(mi,n2) = ham — kr) f (1, 2) 
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For every value 7 of the variable n;, we have 
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Thus, obtaining each output y in row nz requires 2 addition 
and one division by the integer constant m x n, except for 
the first y of the row corresponding to n; = O, m additions 
and a division by an integer constant are required. On the 
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Figure 6: Computation of y(z + 1,2) from y(2,n2). 


other hand, our VFMF technique requires a memory variable 
to store S; obtained at each step in order to be used in the 
next step. Figure [6] illustrates the recursive expressions of y’s. 

To obtain entirely the output signal y based on the input 
image signal x, our VFMF technique needs [(2 x M x (N+n— 
1))+M xnj]+[2x (M+m-—-1)x(N+n—-1)+(N+n—-1) xm] 
additions and (M+m-—1) x (N+n-—1) divisions by integer. 
Since the size of the original image is significantly larger than 
the size of the filter G.e., M,N >> m,n), then approximately 
(4 x M x N) additions and M x N divisions by integer are 
required. 


V. EXPERIMENTAL RESULTS 


We test our VFMF technique on images of various sizes 
using many filter sizes in Microsoft Visual C++ MFC software 
and compare our results with those for the standard, and 
separation implementations. Figure [7]8] and |9] presents the 
results of our experiments. 
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Figure 7 


Figure |7| shows that the SSC technique and our VFMF 
technique provide a huge improvement in speed as the size 
of the filter grows. In addition, it also shows that the standard 
implementation requires a significant amount of computation 
time as the size of the filter grows, which easily reaches 
an unacceptable threshold for applications (typically real-time 
applications) that require high response rates. 
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Figure 8 


Figure[8|]shows the large differences in the computation time 
taken he SSC technique and our VFMF technique as the size 
of the filter grows. 
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Figure 9] shows amount of time to perform our VFMF 
technique of input images of various sizes using different filter 
sizes. We note that the amount of time required to implement 
our VFMF technique is almost invariant under the change in 
size of the filter. 


VI. CONCLUSIONS 


We propose a new very fast mean filtering (VFMF) tech- 
nique for 2-D mean filtering. To perform a 2-D mean filtering 
of an Wx N image with an m x n mean filter, our VFMF 
technique needs approximately 4 x M x N additions and 
MN divisions by integer. The gain in the performance of our 
VEMF technique is significant, since the major additions and 
divisions required by our VFMF technique depend strongly 
on the size of the input image but not on the mean filter size. 
Our experimental results on images of various sizes using a 
mean filter with different sizes demonstrate that our VFMF 
technique is consistently faster than the other techniques, 
especially when the sizes of the original image and the filter 
are large. 

The application of our VFMF technique in the form of 
software can help increase the processing speed, bringing 
high feasibility for software applications requiring real-time 
processing. The application of our VFEMF technique in the 
form of hardware, i.e., a processing circuit can help reduce 
the calculation elements (nodes), thereby positively affecting 
the design and production costs. 

Further, our VFMF technique can be expanded for multi- 
dimensional signal filtering problems by spatial averaging 
filters, which cannot be done by techniques and due to 
their simple and limited approach in two dimensions. 
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